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Radiation by the atoms of a resonant medium is a cooperative process in which the medium par- 
ticipates as a whole. In two previous papers [1, 2], we treated this problem for the case of a medium 
having slab geometry, which, under plane wave excitation, supports coherent waves that propagate 
in one dimension. We extend the treatment here to the three-dimensional problem, focusing princi- 
pally on the case of spherical geometry. By regarding the radiation field as a superposition of electric 
and magnetic multipole fields of different orders, we express it in terms of suitably defined scalar 
fields. The latter fields possess a sequence of exponentially decaying eigenmodes corresponding to 
each multipole order. We consider several examples of spherically symmetric initial excitations of a 
sphere. Small uniformly excited spheres, we find, tend to radiate superradiantly, while the radiation 
from a large sphere with an initially excited inner core exhibits temporal oscillations that result 
from the participation of a large number of coherently excited amplitudes in different modes. The 
frequency spectrum of the emitted radiation possesses a rich structure, including a frequency gap 
for large spheres and sharply defined and closely spaced peaks caused by the small frequency shifts 
and even smaller decay rates characteristic of the majority of eigenmodes. 

PACS numbers: 42.50.Nn,42.50.Gy,42.25.Kb,42.25.Gy 



I. INTRODUCTION 

A quantum of light emitted within a finite medium 
made of identical atoms cannot easily escape. It may 
suffer a large and undeterminable number of coherent 
absorption and reemission processes before reaching the 
boundary of the medium, and even there is subject to the 
hazard of internal reflection. The emission of a quantum 
by a single atom within such a medium, in other words, is 
inevitably a process in which the entire medium partakes 
coherently. When more than one atom is excited initially, 
they tend to radiate cooperatively and that sort of coher- 
ent emission process is often referred to as superradiance 
[3-8] . The time dependence of the energy emitted by the 
medium in all these cases may differ considerably from 
the exponential decay characteristic of the isolated atom, 
and its frequency spectrum may differ substantially from 
the familiar Lorentzian form. In the present paper we 
address ourselves to the analysis of these effects for the 
important case in which the medium is spherical in shape 
and all the inhomogeneities that destroy coherence are 
assumed negligible. 

The medium we envisage consists of identical atoms, 
all with the same electric dipole resonance at a (renor- 
malized) frequency ujq. Since they are never more than 
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weakly excited, the atoms may be replaced, in effect, 
by harmonic oscillators of the same frequency. We as- 
sume them to be distributed densely enough that many 
are present in each cubic wavelength (27rc/cLio)'^, and 
smoothly enough to permit treating the medium as a 
continuum. We have called this idealized model of a reso- 
nant and isotropically polarizable medium polarium, and 
have discussed a number of its behaviors in two previous 
papers, I [1] and II [2]. In those papers, we discussed 
the emission and propagation of radiation in the essen- 
tially one-dimensional context of parallel slab geometry. 
The plane waves that are radiated, we showed, can be 
regarded as a superposition of contributions from a se- 
quence of mutually orthogonal polarization modes that 
decay exponentially and have easily calculated proper- 
ties. 

The modal decomposition of the polarization in the 
one-dimensional medium revealed in I a complex struc- 
ture for the time dependence of its decay. An oscilla- 
tory exchange of energy takes place, in effect, between 
the coherently coupled radiation and polarization fields. 
The radiated spectrum consequently exhibits an elabo- 
rate structure of narrow peaks and dips corresponding 
to mutually interefering resonant amplitudes contributed 
by the various modes. A gap is also present in the spec- 
trum, corresponding to a band of frequencies in which the 
implicit dispersion law of the medium suppresses prop- 
agation. Not surprisingly, these features are also found 
to play an important role in the discussion of reflection 
and transmission of an externally incident plane wave 
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by a slab shaped medium that we undertook in II. By 
changing the angle of incidence and the polarization di- 
rection of the incident wave, we could tune and alter the 
spectral dependences of the reflection and transmission 
coefEcients in predictable ways. 

One of our early tasks in analyzing the problem in 
spherical geometry will be to find the appropriate set 
of exponentially decaying polarization modes. These 
must obey conditions that assure the transverse char- 
acter of the radiated fields. In I, transversality was 
easily secured by dealing only with fields uniform over 
planes perpendicular to the axis of propagation. In 
II, where plane waves could be obliquely incident upon 
the slab shaped medium, a more careful treatment of 
the separated Cartesian field components was required 
since induced surface charges and currents complicate the 
boundary conditions. 

For the spherical geometry we find it particularly con- 
venient to introduce two familiar angular momentum pro- 
jection operators to separate the full vector field problem 
into its electric and magnetic multipole parts. Each of 
these parts can be further resolved into a succession of 
spherical harmonic components. The magnetic multipole 
fields separated in this way automatically obey the re- 
quired transversality conditions. The electric multipole 
fields, on the other hand, still require some considera- 
tion of surface charges and currents in order to secure 
transversality. The most convenient feature of the an- 
gular momentum decomposition is that it involves only 
scalar functions which are in effect the radial components 
of the electric and magnetic fields. We may express these 
components in terms of mode functions that decay expo- 
nentially with time in each spherical harmonic order. 

We begin in Sec. II with the elementary example of 
radiation from a small uniformly excited sphere. This 
problem, which can be solved directly without the use of 
angular momentum operators, provides physical insights 
into the coherent radiation problem, which are useful in 
treating the problem of larger spheres. In Sec. Ill, we 
present the equations of motion for the radiation field and 
polarization. The complications of maintaining a trans- 
verse displacement field everywhere are handled here by 
means of the two angular momentum projection opera- 
tors. Their use leads to equations of motion for the elec- 
tric and magnetic fields and for their polarization sources. 
The magnetic multipole radiation, which is analytically 
simpler to treat, is discussed in Sec. IV, while a full 
treatment of the electric multipole radiation is begun in 
Sec. V. The magnetic multipole radiation field obeys 
simple outgoing- wave boundary conditions in each angu- 
lar momentum order. We derive the eigenvalue equation 
that results from these boundary conditions and briefly 
discuss certain important properties of the exponentially 
decaying eigenmodes. The electric multipole radiation 
problem requires somewhat more involved outgoing- wave 
boundary conditions that can be most simply treated 
by considering the differential equations obeyed by the 
electric and displacement fields inside and outside the 



medium. When the medium excitation is restricted to 
initial polarizations that are oriented along a fixed direc- 
tion and have a radially symmetric but otherwise arbi- 
trary amplitude, the spherical medium, regardless of its 
radius, radiates as a pure electric dipole, as we show in 
Sec. V. In the same section, the associated electric-dipole 
eigenvalue problem is developed in the more general con- 
text of electric multipoles of arbitrary order, and the 
eigenvalue equation is analytically solved for the eigen- 
values in several important limiting cases. In Sec. VI, we 
establish the orthogonality of electric multipole modes of 
an arbitrary order. In Sec. VII we discuss the temporal 
and spectral characteristics of the electric dipole radia- 
tion from an interesting example of a spherical system, 
one in which a uniformly excited spherical core is sur- 
rounded by an initially unexcited spherical shell. Finally, 
in Sec. VIII, some concluding remarks about the prob- 
lem of coherent transport of resonant excitations treated 
here axe presented. 

II. AN ELEMENTARY EXAMPLE: 
RADIATION FROM A SMALL UNIFORMLY 
EXCITED SPHERE 

It will be useful to begin our analysis by discussing 
the radiation by a spherical medium of radius R much 
smaller than the reduced wavelength l/feg = c/ojq, i.e., 
/3 = koR « 1. The problem is simple enough to afford 
elementary access. It will furnish a valuable example for 
later reference. 

The atoms of the medium we call polarium are assumed 
to be distributed with a uniform density no- The tran- 
sition matrix elements of their electric dipole moment 
vectors fl are assumed to be randomly oriented so that 
the medium is isotropic, and its inducxxl electric polar- 
ization is always parallel to the inducing field. Then, as 
we have shown in deriving Eq. (9) of I, the positive fre- 
quency part of the polarization field P^+^(r, for such 
a medium is driven by the positive frequency part of the 
electric field (r, t) through the relation 

(|+i.o)p^-)(f,*) = ^^^(-)(f,*). (1) 

In the absence of the electric field, the polarization varies 
in time at any point as exp(— iwof). 

By expressing and p(-l-) jn terms of their slowly- 
varying envelopes £ and V: 

(2) 

we may rewrite Eq. (1) as 

We shall assume that the polarization is spatially uni- 
form and given by exp(— iwoOi where £ is a unit 
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vector along the z— axis. In that case, the fields gener- 
ated by the uniform polarization, although rapidly os- 
cillating in time, have the familiar electrostatic spatial 
dependence within the near field zone. Thus the electric 
field within the sphere is uniform and parallel to the po- 
larization. The magnetic field within the medium is of 
relative order koR « 1, and is thus negligible in this, 
the long wavelength limit. 

Such a uniformly polarized small sphere decays super- 
radiantly, as we shall see presently. If Aq is the exponen- 
tial decay constant for this mode, then the electric field 
within the medium, which we denote by , and the po- 
larization Poexp(— Ao^) are formally related, according 
to Eq. (3), by 
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Outside the medium, the electric field E-' is that of a 
point dipole po, equal in value to the dipole moment of 
the sphere. 
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and located at its center. In the near field zone, the 
electric field thus assumes the familiar elecrostatic spatial 
dependence 



-> _ 3{po ■ r)f 



E> = 
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(6) 



The eigenvalue Aq is determined by requiring that in 
the long-wavelength limit the electric field in the interior 
of a uniformly polarized small sphere be minus a third of 
its polarization, E^ = — (1/3)P. In view of Eq. (4), this 
requirement yields the following expression for Aq: 
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Because the eigenvalue (7) is purely imaginary, it rep- 
resents merely a frequency shift for the field and polar- 
ization. The sphere does not radiate in the limit /3 0. 
For small but finite /3, the sphere does, in fact, radiate. 
To calculate the decay rate ReAo, we resort to an anal- 
ysis based on the rate at which a uniformly polarized 
sphere, with its polarization oscillating at frequency wq, 
radiates energy. A small sphere of this kind radiates as 
an oscillating point dipole [9] at the time-averaged rate 
[10] 
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Equivalcntly, the rate at which energy is lost from the 
dipoles of the spherical medium must be equal to the 
rate at which work is done by them on the field. The 
rate of work done by the ith dipole of dipole moment pi, 
when averaged over the fundamental period of oscillation. 



is just — 2Re ~p ■ E, where fi^ is the time derivative of 
the dipole moment. Therefore the rate at which work is 

done by the sphere is -2Re "P ■ E{ATrR^/3). Since the 
polarization and field both oscillate at a frequency close 
to LOo, the preceding expression is essentially the same as 



W=^^iJoImP*-E. 



(9) 



When the relation (1) between the field and polarization 
is used to eliminate E from Eq. (9) , the resulting expres- 
sion for W is 

Re Aq I ,^i2 



(10) 



An explicit expression for the decay rate Re Ao is now ob- 
tained by equating Eqs. (8)) and (10) and using Eq. (5), 



ReAo 



2 nol/Ip 



(11) 
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The rate (11) at which a small uniformly polarized 

sphere decays radiatively may be expressed as the prod- 
uct of the number of atoms, N = no(47ri?'^/3), and 
the Wigner-Weisskopf intrinsic decay rate of each atom. 



-1 = im^k^Jimnh): 



Re Ao = Nt- 



(12) 



Each dipole in an assembly of N identically prepared 
and coherently coupled atomic dipoles emits radiation at 
a rate that is N times the rate with which it would spon- 
taneously radiate when isolated from the others. Such 
enhanced decay rates are characteristic of superradiant 
emission, a process that is often considered in the context 
of much stronger excitation of the radiating medium, e.g., 
when all of the atoms are fully excited in the initial state. 
In these more general examples the emission process can 
only be described adequately by means of nonlinear equa- 
tions. The present problem, by contrast, is considerably 
simpler due to its linearity in the fields and polarization. 
Because of the coherent initial preparation of the atomic 
dipoles, the essential coherence always remains present 
in the emission process. 

It is worth recalling here that Hartmann and collabora- 
tors [11, 12] have made an important criticism of Dicke's 
elementary theory of superradiance in many-atom sys- 
tems. They have pointed out that the electric dipole mo- 
ments induced in different atoms will interact strongly 
via the familiar dipolc-dipole interactions and lead to 
spatially dependent shifts of atomic energy levels. These 
differing level shifts can bring about relative dephasing of 
difii'erent parts of the oscillating polarization distribution 
and thus some; breakdown of the cooperative character of 
superradiant emission. We see no evidence of this sup- 
pression in the radiative rate given by Eq. (12). Indeed, 
the way in which we have treated the interaction of each 
atom with the field implicitly includes the effects of all 
dipole-dipole interactions. Their total effect does not in- 
hibit superradiance, at least for the linear problem of 
radiation from a small uniformly polarized sphere. 
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III. FORMULATION OF THE GENERAL 
PROBLEM 



We now turn to the general problem of radiation by 
an arbitrary excitation of a spherical medium of arbitrary 
radius. The resonant interaction of the; polarium medium 
with the electromagnetic field is described by Eq. (1) and 
the Maxwell wave equation 

-Vx(Vx£H))-^*fiH» = ^^pH). (13) 

Because of the identity 

-V X (V X E^+^) = V^^(+) - V(V • ^(+)), (14) 

Eq. (14) contains an explicit gradient term, V(V • E^~^^), 
which enforces the transversality of the total displace- 
ment field, D = E + P. We shall see later that this term 
greatly influences the character of the radiation and most 
particularly when the spheres are small compared to the 
wavelength of radiation. 

The spherical geometry of the radiation problem is best 
approached by decomposing the radiation field into its 
electric multipole (EM) and magnetic multipole (MM) 
components, and introducing appropriate scalar func- 
tions to describe them. A simple way to exhibit this de- 
composition without introducing the full panoply of vec- 
tor spherical harmonics is to use two operators [8. 9] that 
arc simply related to the quantum-mechanical angular- 
momentum operator, L = —if x V. Let us consider the 
action of the operators L- and L • Vx on Eqs. (1) and 
(13). Because both these operators annihilate the gradi- 
ent term inside the double curl when identity (14) is used 
and because they commute with the Laplacian, Eq. (13) 
simplifies to an inhomogeneous scalar wave equation of 
the general form 

where the symbols, <p and r], denote the functions 

(/> = L- P(+) and = i- £;(+)) (16) 
or alternatively 

= V X and ?7 = L- V X ^W). (17) 

That L ■ E^+l and L • V x describe the radial 

components of the magnetic and electric fields and thus 
the MM and EM fields, respectively, is immediately evi- 
dent when the vector triple products are rearranged, and 
the Faraday and Maxwell- Ampere laws are introduced as 



follows: 

L ■ ~(f X V) • ^(+) = r • (V X E(+)) 
dt 

L-W X ~L • ~ f • (V X SW) 

r^r-^^EM r^yf.E^+), (ig) 

The final step in each of the two relations in Eq. (18) 
has employed the assumption of quasi-monochromaticity, 
which permits replacing time differentiations of the 
positive-frequency parts of the electromagnetic field by 
the multiplicative factor —iujQ. The validity of this as- 
sumption is assured by the resonant character of the ra- 
diative interactions in the medium. 

By contrast with Eq. (13), Eq. (1) is formally un- 
changed when the spatial operator L- is applied, 

(I +ia;o) mt) = '-^Vir,t), (19) 

where <j) and rj are given by Eq. (16). Under the L ■ Vx 
operation, however, the sharp drop-off at the surface of 
the otherwise uniform medium density contributes to the 
right-hand side of Eq. (19) a surface singularity [13], 
which is generated by oscillating surface charges. Such 
surface polarization charges and currents must be present 
for EM radiation for which both the electric field and po- 
larization have nonvanishing radial components. A self- 
consistent approach that treats these surface singularities 
correctly but implicitly is based on matching on the sur- 
face the appropriate components of the electromagnetic 
field inside the medium to those outside. 

The L- operation commutes with the density which has 
only a radial step-function form. Eq. (19) is thus valid 
both in the interior and on the surface of the medium. 
That is to say, for MM radiation, the electric field and 
polarization are both purely transverse, and no surface 
charges are present. 

For the problem of radiation from an initially excited 
medium with no externally incident fields, Eq. (18) ad- 
mits the familiar sort of retarded integral solution for 77: 

= - 4 / : 1^ ^ dr'. (20) 

c2 y 47r|r-r'| 

By expressing -q and (j) in terms of their slowly- varying 
envelopes £ and V: 

r,{f, t) = Sir, t)e-'^°\ ^{r, t) = V{r, t)e-'^-\ (21) 

and dropping the time derivatives of V , which may be 
assumed small, we may reduce Eq. (20) to the form 
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Because Eq. (19) fails to include radiating surface cur- 
rents, the description provided by Eqs. (19)) and (22) is 
not complete for EM radiation. But these equations do, 
however, describe properly MM radiation which has no 
surface sources for a spherical radiator. 



IV. MAGNETIC MULTIPOLE RADIATION 
AND THE ASSOCIATED EIGENVALUE 
PROBLEM 

The use of the envelopes (21) in Eq. (19) leads to the 
equation of motion for the polarization multipoles. 



(23) 



By eliminating the field multipoles £ between Eqs. (22) 
and (23), we secure the integral equation for the polar- 
ization multipole fields 

d^,^,. jMlW g /•e^'^ok-.-l 



47r|r — r 'I 



df'. 



(24) 

The integrand on the right side of Eq. (24) requires 

evaluating the envelope function V at the retarded time 
t — \ f ~ r'\/c, but if its temporal variation is sufficiently 
slow it remains accurate to neglect the retardation in V 
and write 

-P(f, t) = j ^-^^V if', t) df'. (25) 

The more important effects of retardation are still re- 
tained in the exponential function in the integrand. This 
approach, we have called the rapid-transit approximation 
in I, assumes only that the slowly varying amplitude V 
does not change appreciably during the passage time of 
a wave through the medium. The approximation is not 
essential to our treatment of the problem, but it greatly 
simplifies the analysis, and we shall therefore employ it 
in exploring the behavior of the system. 

Let us consider an expansion of V{f,t) in spherical 
harmonics, defined according to the convention employed 
in Ref. [13], 



V{r,t) = J2Pem{r,t)Yemi^). 



(26) 



Substitution of this form into Eq. (25), followed by a use 

of the identity 

'^[^37^ = ^koY,Je{kor<)h''l\kor>)Yer,,{n)Y;^{n% 
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and integration over the solid angles f1' of the vector 
r ', together with a use of the orthonormality and linear 
independence of the various spherical harmonics yields 
the result 

xje{kor<)h^P{kor>y^dr'. (28) 

Here je and h^P are spherical Bessel and Hankel functions 
of the first kind, and r<(r-*) is the smaller (larger) of 
the two radial distances r, r'. That each multipole order 
separates from all others is a consequence of the spherical 
geometry. All multipoles that are initially unexcited thus 
remain unexcited at all later times. 

We now look for solutions of Eq. (28) that have a purely 
exponential time dependence. 



-At 



(29) 



For such solutions, Eq. (28) reduces to the homogeneous 
Predholm integral equation: 

no|/?|2/c3 rR 



X jt{kQr<)hP [kor^y^dr' . 

(30) 

This equation contains the radiative boundary condition 
at the spherical surface r = R 



(31) 



(27) 



together with the condition that PemX remains finite at 
r = 0. These conditions restrict the form of the po- 
larization function Vemxir) and the values of the decay 
constant A. 

The possible values of A, the eigenvalues, form a dis- 
crete, infinite set of complex numbers. By invoking the 
symmetry of the kernel of the integral equation (30) un- 
der the interchange r -f-)- r', we can establish, as in I, both 
the positivity of the real part of each eigenvalue and the 
orthogonality of the eigenfunctions, 

l-R 

/ remx{r)rimX'{rydr^6xx'- (32) 
Jo 

Equation (32) can be used to define a normalization inte- 
gral and to secure an orthonormal set of eigenfunctions. 
Further properties of the eigenvalues and eigenfunctions 
follow from the bilinear expansion of the kernel in terms 
of the orthonormal eignfunctions, as we showed in the 
context of the one-dimensional problem [1]. 

To solve explicitly for the eigenvalues A and eigen- 
functions Vimxif), we employ the fact that the kernel 
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ji{kQr^)h^/\kor'') is a Green's function, 

2 d 



\ dr"^ r dr 



+ 



'^0 „2 



]dk^r<)h^l\kor>) 



(33) 

This enables us to convert Eq. (30) into a differential 
form, 



d? 2d_ 

(ly.2 J. fly. 



Vemx{r) = 0, (34) 



where 



7a = ko\ 1 



mo|/?P/(3/t) 
A 



(35) 



may be regarded as the propagation constant for the 
mode inside the medium. 

The most general solution of this equation that remains 
finite at r = has the form 



(36) 



The radiative boundary condition (31) at r = R is 
equivalent to the following equality involving logarithmic 
derivatives: 



h^Pik^R) 



(37) 



which immediately leads to the eigenvalue equation 

Ixj'tilxR) _ kahf^\koR) 



m{i\R) 



hf\kQR) 



(38) 



where the prime superscript denotes first derivatives of 
the respective functions with respect to their arguments. 



V. EXCITATIONS OF SPHERICALLY 
SYMMETRIC AMPLITUDE AND ELECTRIC 
DIPOLE RADIATION 

If the initial polarization present in the medium has 
complete spherical symmetry, it must point everywhere 
in the radial direction with an amplitude that has no an- 
gular dependence. Such polarizations, being purely longi- 
tudinal, cannot radiate at all. The radiation of transverse 
waves requires that the spherical symmetry be broken. 

We consider radiation by initial polarizations of the 
spherical medium that have a uniform direction through- 
out the sphere but spherically symmetric amplitudes. If 
we take this direction to be the z axis identified by the 
unit vector z, we may write 



p(+)(f*,t = 0) = zp{r). 



Since z commutes with L which annihilates any spher- 
ically symmetric function, the scalar product L ■ p(+) 
vanishes initially and, according to Eq. (25), at all subse- 
quent times as well. Because of Eq. (23), L ■ there- 
fore also vanishes at all times. The initial excitation (39) 
thus cannot radiate MM fields. 

The operation of L • Vx on Eq. (39), on the other 
hand, produces a nontrivial result. Because of the den- 
sity discontinuity at the boundary, this operation involves 
a surface singularity, as we noted earlier. It leads, fur- 
thermore, to a finite result within the medium, that may 
be shown after simple algebra to involve the derivative of 
the amplitude p(r), 



L - V X p(+)(f,0) 




(40) 



(39) 



An excitation of the form (39) thus radiates only an elec- 
tric dipole contribution of the {t = l,m = 0) order. All 
other multipoles remain unexcited. 

It is interesting to compare our coherently excited ini- 
tial state with an entangled single-excitation state of an 
extended medium of identical two-level atoms envisioned 
by Scully and collaborators [14, 15] as either having been 
prepared by a swept-wave excitation or being intially in 
the symmetric Dicke state. While the entangled atomic 
quantum states differ from our coherent initial excita- 
tion in essential ways, they all share the important at- 
tribute of initial coherence, as easily confirmed by the 
nonzero off-diagonal matrix elements of the density op- 
erator for the single-excitation state. We claim that it 
is this spatially extended initial coherence, not entangle- 
ment per se, that is fundamentally responsible for cooper- 
ative radiation processes such as superradiance and sub- 
radiance. The absence of sub-radiant emission for the 
swept-excitation state is a result of a coherent phasing 
of the emitted photon that yields an enhanced emission 
rate in the forward direction. Both for our problem and 
the symmetric Dicke state, however, modes in which the 
atoms cooperate to trap coherent excitation and release 
it only weakly are also possible when /3 >> 1. But unlike 
the single-excitation state analysis, which is based on a 
scalar- field approximation [14, 16], our exact vector-field 
treatment accounts fully for the polarization and angular 
distribution of the emitted radiation. 



A. Solution Procedure 

We first develop a general approach to the electric mul- 
tipole radiation problcini. and then apply it to the; special 
case of elctric dipole radiation we have just discussed. 
Instead of formulating the problem in terms of an inte- 
gral equation, as we did for the MM radiation, we begin 
with the differential equations that describe the radia- 
tion fields inside and outside the spherical medium. In 
the rapid-transit approximation, which we use, the retar- 
dation of the slowly varying amplitudes is negligible, and 
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that reduces Eq. (22) to the simple form 

iko\f-f'\ 



Sif,t) 



rdr'. 



(41) 



The integral equation (41) can be turned into a differ- 
ential equation by operating on both sides of it with the 
operator (V^ + ^q), 

(V^ + kl)£{r, t) = -klV{r, t). (42) 

Wc look for solutions of Eqs. (23) and (42) within the 
medium that are expansions in spherical harmonics of 
the form 

P<(f,t)=^P<„(r,t)F,„(i1), 



S<{r,t)=Y^£Ur,t)YU^). 



(43) 



The superscript < on the variables V and £ is here used 
to identify the polarization and field inside the medium. 
On writing the Laplacian as 



d 



2 d 



Qj.2 ^ J- Qf 7-2 



(44) 



and noting that the spherical harmonics ¥(„. are mutu- 
ally orthogonal eigcnfunctions of the operator with 
eigenvalues £{i + 1), Eq. (42) separates into individual 
equations, one for each spherical harmonic, 



2d_ 

r dr 



«,Q 



= -klVUr,t). (45) 



The evolution of the polarization within the medium, de- 
scribed by Eq. (23), also separates similarly. 



d 



(46) 



Outside the medium, i.e., for r > R, the polarization 
vanishes identically, and the EM field obeys the free- 
space wave equation obtained by setting equal to 
in Eq. (45). Denoting the exterior field by the super- 
script >, we may thus write for r > R 



(d^ 2d_ 

\ dr'^ r dr 



0. (47) 



Our radiation problem thus separates both inside and 
outside the medium into radiation by the individual mul- 
tipoles. 

The interior field and polarization are coupled to the 

exterior field by the continuity of the normal component 
of the displacement field Df^m and of the tangential com- 
ponents of the electric field at the boundary, r = R, 



and 



(48) 



(49) 



B. The Eigenvalue Problem 

Because of the separability of the space and time vari- 
ables in Eqs. (45)- (47), they admit exponentially time- 
dependent solutions that maintain their shape and are of 
the form 



-xt 



(50) 



For such solutions, Eq. (44) defines the constitutive rela- 
tion between the polarization and the field 



'PemX = -i- 



3hX 



-£Lxir)- 



(51) 



When combined with Eq. (45), this relation leads to the 
equation for £^x- 

92 



2d_ 

r dr 



+ 



^(€+1) 



^Lxir) 



0. 



A general solution of Eq. (52) that is also finite at r 
must, as we have seen earlier, take the form 

^Lxir) 



(52) 
= 



Auilxr) (53) 

within the medium. Outside the medium, the fields con- 
sist of purely radiative, out-going waves. They are de- 
scribed by a solution of Eq. (47) of form 



^f^x{r) 



Bhf\kar) 



(54) 



From these field forms, we obtain, as we now show, the 
electric, magnetic, and displacement fields throughout 
space. 

Returning to the definition of rj in Eq. (17), of which 
£ is the slowly varying envelope, we note that the full 

electric field Efm\ and the radial component £e,mxyim of 
its envelope are related as follows: 

L • V X ^,„A - £trnx(r)Ytm{^)e-^'e-''^-K (55) 

From Faraday's law, which, in the slowly-varying en- 
velope approximation, takes the form, V x Ei^X ~ 
ikoBimX, we have 



L ■ BimX = -rr^lmx{r)Ytm.{^y' 

iko 



-Xt^—iujot 



an equation that has the simple solution [9 
1 



Bemx = 



ikol{l + 1) 



£cmx{r)LYem{^)e-'''e 



Xt —iujQt 



(56) 



consistent with the requirement that B has no radial 
component for pure EM radiation. If we now use 
the Maxwell-Ampere law in which the time derivative 

(dD/dt) is approximately set equal to —ioJoD, we obtain 
the corresponding displacement vector 

DlmX ~ V X B^mX 

iko 



k^l 



■ + 1) 



V X £imx{r)LYem{^) 



-{X+iuio)t 



(57) 
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Equation (57) also describes the electric field in the free 
space outside the medium. Finally, the electric field E^^^ 
within the medium is obtained by means of Eq. (51), 
which yields 



D 



imX 



3hX 



imX- 



This equation is easily solved for Ef^^, 



^2^lmX^ 
IX 



(58) 



an expression in which use has been made of the defini- 
tion (35) for 7a. The ratio {'y\/ko)'^, the electrical per- 
mittivity of the medium for the propagation of the {imX) 
mode, naturally connects the electric and displacement 
fields within the medium. 

We are now in a position to apply our boundary condi- 
tions (48) and (49). Taking the scalar product of Eq. (57) 
with r and noting that 

f-{VxF) = {fxV)-F = iL-F 

for an arbitrary vector field F and that L commutes with 
any function of the radial coordinate r, we find the rela- 
tion 



r-DimX = Ta^t-mxir) 



£{£ + 1) 



-r2^enix{r)Yirn{^)- 

(59) 

In view of this relation, the condition (48) amounts sim- 
ply to the continuity, at the boundary r = i?, of the 
electric field amplitude £gmX given by Eqs. (53) and (54), 



(60) 



The expression for the tangential components of EgmX 
is more involved. By taking the vector product of 
Eq. (57) with f and using the rule A x {B x C) = 
{A ■ C)B — {A ■ B)C while allowing for the correct spatial 
derivatives to be taken, we have 



r X DimX = 



k^i{i + 1) 



riV£emx{r)LiYer. 
d 

r -g;^£emx{r)LYim 



in which it is understood that the repeated vector- 
component index i is summed over its three values. Not- 
ing now that fVg = V{fg) — gVf and that rjL, = f-L = 
0, we have the expression for f x D^mX in terms of the 
amplitude £e,mX 



r X Demx = - 



f,2g(l^ ^-^ i^emx + rS'emxir)] LYem- (61) 



But since E>^^ 



Dfmx and from Eq. (58) E<^^ 



iDf^^, the boundary condition (49) may be written 



as the matching condition 

T2 [^fmX + ^fmx(^)] = i^lmX + ^i'mxi^)] ' (^2) 

When Eqs. (53) and (54) are substituted into this rela- 
tion, we secure the second condition that our solutions 
must obey. 



— MjxR) + IxRj'tilxR)] 

IX 



B 



h^l\koR) + kQRh\"''{kQR) . (63) 



The two conditions (60) and (63) must be met simulta- 
neously, and we have, by taking their ratio, the eigenvalue 
equation 



where 



/3 = fco-R and x = ^xR = P\ 1 — ^ 



3fiA 



(64) 



(65) 



This eigenvalue relation applies to radiation from an elec- 
tric multipole of arbitrary order {£,m). For our special 
initial condition (39), we need, however, only consider 
the azimuthally symmetric, electric dipole radiation with 
^= l,m = 0. 



1. Electric Dipole Radiation 



Since 



{xji{x)y = sina; -h 



cos X sm X 



X x^ 



and 



(M'^(^))'=(-^ + ;^ + ^ie-' 



(66) 



Eq. (64) reduces to the form 

sin a; 1 i i 

sin a; — a; cos a; a;^ (3 (3{(3'^ + i(3)' 

from which, by a simple transposition and division of 

both sides of the equation by sin a;, the following tran- 
scendental form for the eigenvalue condition results: 



X cot X = 1 



/3a;2 



(67) 



Although still complicated in appearance, the form (67) 
for the eigenvalue equation permits an asymptotic anal- 
ysis of the roots x = jxR- and hence of the eigenvalues 
A - in the limits of large and small radius, /3 >> 1 and 
/3 « 1. 
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2. Electric Dipole Radiation from Small Spheres, /3 << 1 

The superradiant mode that we discussed in Sec. II is 

but one of a sequence of exponentially decaying radiation 
modes that one can discuss for a small sphere. As we shall 
see presently, all of the other modes, however, radiate 
quite weakly by comparison. 

The superradiant mode has a propagation constant 
X = jxR that is of order f3, as we shall see shortly. For 
/3 << 1, we may expand the left hand side of the eigen- 
value equation (67) in powers of (3, retaining only terms 
of the most significant order separately in the real and 
imaginary parts on the left hand side to obtain the ap- 
proximate equation 



X cot X : 



— t 



(l-a;V/32) (l-a;V/32)2- 



If we now expand the left hand side in powers of x and 

keep only the two most significant terms, we may further 
approximate the eigenvalue relation for small |a;| as 



-2r + i 



(l-a;V/32 



an equation that may be solved by an iterative procedure. 
To the lowest significant order in its real and imaginary 
parts, the solution assumes the expression 



-2/32 - 2i(3^ 



(68) 



By employing the relation (65) between x and the eigen- 
value A, we then obtain 



A; 



27h 



(69) 



This expression for the complex decay constant of the 
fundamental mode coincides with that obtained in Sec. II 
by means of physical considerations. 

For all other eigenmodes, x is of order 1 or larger, so 
we can neglect the term /3 in comparison with ix"^ in 
the denominator of the right hand side of the eigenvalue 
equation (67) and obtain the approximate expression. 



X cot a; = 1 



'l-il3-l3'^' 



By expanding the denominator for small /3 and keeping 
only terms of the most significant order separately in the 
real and imaginary parts of the resulting right hand side 
of the preceding equation, we may approximate it as 



a; cot a; = 1 -|- ip^. 
Solutions of Eq. (70) of form 

n - - j IT + e„, n 



1,2, 



(70) 



(71) 



where |e„| << 1, may be found, as we see by substitut- 
ing Eq. (71) into Eq. (70). To the lowest order in e„, 



the following expression for the allowed values of x thus 
results: 



Xn= \ n- 



TT — I 



0'' 



(72) 



From the relation (65) between x and A for any mode, 
we may now solve for the latter. On keeping only terms 
to the lowest order in B separately for the real and imag- 
inary parts, we obtain the nth eigenvalue 



— 



2 nol^l' 



P' 



h {n ■ 
.noiMp 



1/2)4 



3h (n-l/2)27r2' 



1,2, 



(73) 



Both the frequency shift and decay rate for each of 
these modes are suppressed by a factor 0(/3^) relative to 
those for the superradiant mode. The dramatic decrease 

of the decay rates with increasing mode index shows that 
the excitations of these modes remain trapped for long 
periods. 

Even when is of order 1, the behavior of the eigen- 
values, A„, is qualitatively similar to that for /3 << 1 we 
have just considered. We demonstrate this by plotting 
the real and imaginary parts of the first 20 eigenvalues 
for /3 = 1 in Figs. 1(a) and 1(b). The existence of only 
a single superradiant mode is quite clear. The next four 
most strongly decaying modes have decay rates that are 
roughly 50, 500, 2000, and 6000 times smaller. Their 
frequency detunings from atomic resonance, which are of 
the opposite sign relative to that for the single superra- 
diant mode, decrease more modestly, however, so their 
successive differences are quite small. These features of 
the weakly decaying modes are responsible, as we shall 
see later, for a slow oscillatory decay of efficiently trapped 
excitation from a small sphere which is initially excited 
from its center out to only a small fraction of its total 
radius. 



3. Electric Dipole Radiation from Large Spheres, f5 » 1 

When the radius of the sphere greatly exceeds the 

wavelength, >> 1, we may ignore (/3^ + i0)^^ when 
compared to 1, and Eq. (67) simplifies somewhat to the 
expression 



X cot a; = 1 — 



/3a;2 



P + ix'^ 



(74) 



If the inequality |a;p >> P also holds, then the eigenvalue 
equation simplifies still further to the form 



X cot X 1 + iP ip. 



(75) 



This equation is formally identical to the eigenvalue equa- 
tion for the odd modes in the one-dimensional problem 
of resonant propagation inside a slab that we discussed in 
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Decay Rates 




(a) 



Frequency Detuning 
6=1 




10 12 14 16 18 20 

Index, n 



(b) 



FIG. 1. (a) Decay rates and (b) frequency detuning (from 
resonance) for the 20 fastest decaying modes, for /3 = 1, in 
units of nop^/3/(3/t). 



I, provided the slab thickness L is taken to be the sphere 
diameter 2R. This isomorphism between the spherical 
and slab geometries in the limit of large medium exten- 
sions holds only for those modes that have propagation 
constants 7a that in magnitude greatly exceed the cur- 
vature 1/R of the spherical surface, consistent with the 
mathematical requirements \x\^ » P and /3 >> 1. 

Because of the infinitely many branches of the cotan- 
gent function, the eigenvalue equation (74) has infinitely 
many solutions. However, only those solutions for which 
I a; I is comparable to /3 correspond to field eigenmodes 
that, being spatially phase matched to the radiation field, 
are strongly coupled to it, and thus radiate efficiently. 
The eigenvalues of the other, relatively weakly radiating 
modes are analogous to those we have already considered 
for the one-dimensional problem in I, and need no further 
attention. 



When |a;| ~ /3, then \x'^ 



» 



and the right 



hand side of Eq. (74) may be approximated to the lowest 
order in its real and imaginary parts by the expression 



i/3 + {1 ~ 0^ jx^). Now dividing both sides of Eq. (74) 
by a; ~ /3, we may reduce it to the following approximate 
form: 



(76) 



cot a;~z-|-— 1 



where we have only retained the most significant power 

of the small quantity (1 — (3 / t) in the real and imag- 
inary parts separately. By substituting x = —iy and 
cot a: = icothy in Eq. (76), we see that cothy differs lit- 
tle from 1, which implies that y must have a large real 
part and cothy w 1 + 2exp{—2y). With this approxima- 
tion, Eq. (76) can be solved for its multiple roots. They 
may each be labeled by an integer n, with the value 



x„ = -iVn ^{n+ l/4)7r - {i/2) In^. 



(77) 



For consistency with the assumption \x\ ~ /3, we require 
that mr ^ p. The eigenvalues A„ now follow from their 
relation (65) with x„. 



mo I /i IV (3ft) 



l-[(n- 



1/4)2^7/32] 



i[(n-hl/4)7r//32]ln/3 



(/3-n7r)-Fi(l/2)ln/3' 



(78) 



where the last approximate equality follows from drop- 
ping the 1/4 when compared to n, and then setting 
riTT w /3 and (^2 _ ^2^2^ _ - mr). The real and 
imaginary parts of the eigenvalues A„ thus have a simple 
resonant character, similar in form to the frequency de- 
pendence of the imaginary and real parts, respectively, 
of the susceptibility of a resonant dielectric medium. 

In Figs. 2(a) and 2(b), we display for l3 = 1000 the 
real and imaginary parts of the eigenvalues in the reso- 
nant domain, mr 0. These have been obtained both 
by a highly accurate numerical treatment of the exact 
eigenvalue equation (67) and by means of the analyti- 
cal approximation (78). The analytical approximation is 
evidently quite accurate for such a large value of /3. 

For further discussions of EM radiation from large 
spheres we undertake the development of the orthogonal- 
ity properties of the EM eigenmodes. They are somewhat 
less self-evident than the analogous properties of the MM 
eigenmodes we considered in Sec. IV. 



VI. ORTHOGONALITY PROPERTIES OF THE 
ELECTRIC MULTIPOLE MODES 

For the nth mode of the (im) order, which wc shall 

call simply the {imn) mode, the magnetic field B^i^^^ is 

given by Eq. (56), and the electric field E^^^^ essentially 
by its curl, 



S(EM) 
imn 



ikoi{£ + 1) 



£emnir)LYemm, (79) 
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FIG. 2. (a) Decay rates and (b) frequency detuning (from 
resonance) for the 40 fastest decaying modes, for /3 = 1000, 
in units of noH^/3/{3K). 



yields the result 



r<R 



{EM) p{EM) ,^ 



lln Jr=R 



]4_ f ^{EM) ^(EM) 
len Jr<R 



(82) 



Here Faraday's law was employed in the last term to ex- 
press the curl of E^fJ^l^, in terms of the magnetic field 

^e^m/n' °^ {i'm'n') mode. That the last term in 
Eq. (82) vanishes unless i = I' and m = — m' follows 
from the orthogonality of the different vector spherical 
harmonics LY^m in terms of which the magnetic fields of 
the two modes may be expressed by means of Eq. (79). 
By a similar but more detailed argument involving the 
identity [9] 



iV X L = rV^ - V 1 + r 



d_ 

dr 



we may reduce the surface integral in Eq. (82) to a form 
involving J cPfl LYim • LYii^m'i which too vanishes unless 
£ = £' and m = — m'. This proves the orthogonality of 
the vector electric fields for modes with different £ and 
m values. 

A different approach is needed to establish orthogonal- 
ity rules for modes with the same £ value. We begin with 
the field equations that e')'^'' and E^^^J obey, 

-fcoVxBir^+^i^l-f =0, (83) 



{EM) 



(84) 



p{EM) _iko^ ^(EM) 
^Imn ~ -,2 ^ ^ ^imn ■ 
'in 



(80) 



Let us first consider the orthogonality integral for the 
electric fields of two modes p, p' in two different multipole 
orders {£,m) and {£',m'), with £ £': 



E, 



{EM) p{EM) 



■ E^Z'Zdf. 



r<R 



Use of Eq. (80) to replace the electric field E^^^^ in 

terms of the corresponding magnetic field -B^^^^ and the 
identity 



\I ■{AxB)=B-V xA-A-V xB, (81) 

followed by the use of Gauss's theorem to reduce a total 
divergence term to a surface integral over the sphere. 



Taking the scalar product of Eq. (83) with Ef^^, and 
subtracting the result from that obtained on taking the 
scalar product of Eq. (84) with E^^^^ , and integrating 
the difference over the spherical sample we secure the 
result 



{itn Itn') I ^Imn 

lr<R, 



= iko I 

Jr<R 



ff{EM) ^ ^{EM) 



— -tj, 



{EM) 



V xB 



S(EM) 



Im'n' 



dr. 



(85) 



We now use the identity (81) to reexpress each term on 
the right hand side of Eq. (85) in terms of a complete 

divergence, which reduces to a surface integral according 
to Gauss's law, and a curl term, which from Faraday's law 
can be written as a volume integral of B^^'^^ ■ -B^^f^K 
We see in this way that Eq. (85) simplifies to a single 
surface integral, the two volume integrals canceling each 
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other out, 



(7fa Itn'l I ^emn 
lr<R 



Jr=I 



(r X E^EM)y-^iEM) 
if X E^'^^h ■ B^^^^ 



(86) 



The surface integral involves only the tangential compo- 
nents of both the electric and magnetic fields, which are 
all continuous across the spherical boundary. The surface 
integral can, as such, be expressed equally well in terms of 
the same components of the free-space fields in the imme- 
diate exterior of the sphere. However, since in each mul- 
tipole order (£m) the exterior fields are all expressed in 

terms of the outgoing wave solutions h^p {kQr)Yim {^) and 
their derivatives, independent of the mode indices n,n', 
the surface integral on the right hand side of Eq. (86) 
vanishes identically. It follows then from Eq. (86) that if 
lin 7^ IM, then the orthogonality integral vanishes. 



r<R 



(87) 



Combining all of the specific orthogonality relations we 
have discussed so far in this section, we may write down 
the overall orthogonality relation 



e.{EM) ^(EM) 



E. 



r<R 



£mn 



'Sr, 



This orthogonality relation (88) may also be directly 
established by using Eq. (79) and (80) to express its 
left hand side in terms of the mode hmctions Sgmn ^ 
jeiltnT), simplifying the resulting integral by means of 
angular momentum operator identities, and then exploit- 
ing the eigenvalue relation (64) that each mode must 
obey. 

It is worth noting that the requirement m = — m', 

rather than m = m! , for the nonvanishing of the ex- 
pression (88) is a reminder of the symmetric but non- 
Hermitian character of the propagation kernel for the 
electromagnetic field. This nonhermiticity was already 
noted in II in the context of one-dimensional propagation 
when an incident wave with transverse magnetic polar- 
ization was obliquely incident on a slab. 



A. Expansion of an Arbitrary EM Field in the 
Corresponding Modes 

The orthogonality integral (88) makes it possible to 
expand an arbitrary electromagnetic field in a particular 
EM order (L m) in terms of the EM modes, given by 
Eqs. (79) and (80), of the same order. Let V x [f{r)LYtm] 
be an arbitrary field in this order expressed as the mode 
sum 



V X 



f{r)LY,^] = J2 



(89) 



The coefficients /„ are obtained by multiplying both sides 
of Eq. (89) by the mode function V x [j£(7^„'r)ZY^,_m], 
integrating over the spherical sample, and utilizing the 
orthogonality relation (88), 



fn 



{k{lenr),jtiritnr)y 
where the symbol (/, g) defines an inner product, 



(90) 



(/,5) = j dr{y X [/(r)Lr,,_j}-{v x [g{r)LY,m]] ■ 

(91) 

Whenever either / or is a mode amplitude function 
jiilen'r), the inner product (91) can be simplified greatly. 
To see this, we make use of formula (81) in Eq. (91) to 
transform its integrand to a term involving V x (V x 
[g(r)LY£m\), which, when g is a, mode function, is simply 
^f^g{r)L Ye„i; along with a pure divergence term, a term 
that from Gauss's theorem is equal to a surface integral. 
Thus, when g is& mode function, Eq. (91) takes a simpler 
form, 



(/> 9) = Tin / f{r)g{r)LY,,_^ ■ LY.^dr 

Jr<R 

+ j f{r)LYe,_^ x [V x {gLYem)] ■ f d^S. 
Use of the solid-angle integral identity 

d^QLYe^-m ■ LYem = (-1)"+V(£ + 1) 



and the fact [9] that the surface integral in the preceding 
equation is simply {-l)"'+H{£ + l)Rf{R){d/dR)[Rg{R)] 
reduces the inner product (91) to the form 



{f,9) = [RfiR)^ mm 

+lLj^''f{r)g{rydry (92) 

This form of the inner product can also be used as the 
starting point of a simple, direct proof of the orthogonal- 
ity of different modes in a given multipole order. 



VII. RADIATION FROM A LARGE SPHERE 
WITH A UNIFORMLY EXCITED SPHERICAL 
INNER CORE 

Coherent radiation from a large resonant sphere will, 
in general, involve a large number of modes of both mul- 
tipole types and their infinitely many orders. A partic- 
ularly simple situation occurs, however, when an inner 
concentric spherical region of the sphere is uniformly ex- 
cited initially, as shown in Fig. 3. This is a special case 
of a polarization with radially symmetric amplitude and 
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We thus have the result 



FIG. 3. Spherical medium of radius R with an excited con- 
centric core of radius ro = fR, / < 1- 



uniform direction, which emits pure electric dipole radi- 
ation in the (1,0) order. 

Let us assume then an initial polarization of form 



^ ^ I for ro < r < ii. 



(93) 



The electric field that this polarization initially radiates 
has a simple expression in the rapid transit approxima- 
tion [1] in which retardation of slowly varying amplitudes 
is ignored and the second time dcaivatives in Eq. (16) are 
replaced by — w^. In this approximation, the initial elec- 
tric field E<^+\f,0) may be expressed as an integral of 
the scalar product of the tensor propagator 



G(f-r") 



1 _ ^\ piko\r-r'\ 

1 + VV ) ,^ (94) 



and the initial polarization (93) over the spherical sam- 
ple, which reduces to the form 

piko\f-f'\ 

Tl^^df'. (95) 
4n\r-r'\ 



- — —df' X z. (96) 
Anlr-r'l 



The initial field amplitude d{r,0) is the radial function 

that multipUes the spherical harmonic Yiq{Q) in the 
quantity Z • V x E{f, 0). Taking the curl of the left hand 
side of Eq. (95) eliminates the pure gradient term, 

^ r piko\f-f'\ 

V X E{r,0) = fcgPoV / 

J r'<ro 

By using the identity (27) to replace the integrand in 
Eq. (96) by a spherical harmonic sum and then integrat- 
ing it over the angular coordinates, we see that the inte- 
gral in Eq. (96) is a purely radial function, 

pro 

f{r)=iko jo{kor<)h^o\kor>y^dr'. (97) 
Jo 



dff 



V X E(r,0) = kiPoVfir) xz = ^- xz. 

dr r 

Operating on both sides of this equation by noting 
that any function of r commutes with L, and using the 
result 

L • (f X £) = i (f X V) • (r X i) 

= -\r- iz ■ V)f - (f - 5)(V • r)l 

= 2if- z = 2ircosd = 2i\ —rYio{il), 
V " 



we obtain 



L • V X E{f, 0) = d{r, O)Fio(O), 



(98) 



where the initial field amplitude d{r, 0) takes the explicit 
form 



joikor<)hi'\kor>y'dr- 



(99) 

The radial integral can also be performed in closed form 
by means of the indefinite-integral identities [17], 



/ 



joix)x dx = X ji{x), J {x)x dx = X h\ (x), 

(100) 

and its derivative taken by means of the identities 

^io(x) = -ji(x), ±hi'\x) = -h['\x). (101) 

The following final form of the amplitude is obtained in 

this way: 



I Air . 



d(r,0) = 2ktris\ ^Poji{kor<)h^T'\kor>), (102) 

where and are defined to be the smaller and the 
larger of the two radial distances r and ro, respectively. 

The same steps that took us from Eq. (55) to Eq. (57) 
can be employed to write down the displacement field en- 
velope D{f, t) in terms of the amplitude d{r, t) of electric 
dipole radiation {(. = 1), 

'^(^' ~ ^^2 V X [d(r, t)LYw{n)]. (103) 

Inside the medium, the amplitude d{r,t) may 
be expanded in terms of the mode functions 

hilnr) exp(-A„t), 



(104) 



n=l 
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Orthogonality of the mode functions under the inner 
product (92) enables us to solve for the coefficients, 



^ _ (d(r,0),.n(7nr)) 
^ (ii(7nr),ii(7nr))' 



(105) 



The inner products in Eq. (105) contain integrals involv- 
ing bilinear expressions of spherical Bessel and Hankel 
functions of the same order 1, which, as we show in Ap- 
pendix A, can be evaluated in terms of simple trigono- 
metric functions. The inner product in the numerator 
assumes the form, 



(rf(r,0),ji(7„r)) 




while that in the denominator greatly simplifies in the 
limit of large |7n|i?, 

ihilnr),ji{7nr)) « -R, 

so the amplitude rf„ may be expressed as 




jiilnro)- 



(106) 



The electric field envelope E{f, t) has a similar form as 
the electric displacement field D{f,t), 

ZKq 

with the amplitude e(r, t) having a similar mode decom- 
position as d{r, t) , 



e(^,*) = ^npji{'ynr)ex.p{-Xnt). 



(108) 



The mode coefficient e„ for the electric field is obtained 
from dn, the analogous coefficient for the displacement 
field, by dividing the latter by the dielectric constant 
7ra/^o fo'^ '^^^ mode. 



dn 



h2 2 




3\{lnro)- 



(109) 

The resonant character of the excitation of the modes is 
unmistakable in the denominator of the right hand side 
of Eq. (109) only modes with propagation constants 7„ 
that are closely matched to the free space wave vector fco 
are preferentially excited. Also, because ji(7raro) falls off 
rapidly with 7„ for |7„|ro >> 1, only modes with prop- 
agation constants of order 1/ro or smaller in magnitude 
are significantly excited in the preparation of the initial 
core excitation. 



A. Frequency Spectrum of Radiation 

The initial polarization (93) radiates light into the var- 
ious electric-field eigenmodes with amplitudes e„. Apart 



from overall angular and distance dependences, the spec- 
tral amplitude of radiation, a{R,duj), at detuning Suj is 
given by the Fourier transform of the electric field ampli- 
tude e{R,t) at the surface of the sphere, 



a((5w) 



3(i?, t)e'^^'dt = (110) 



n=l 



A„ — iSco 



By substituting for e„ from Eq. (109) and using the rela- 
tion (35) between 7a and A in Eq. (109), we may express 
the spectral amplitude as 



where Ps,^ is R times the propagation constant for a plane 
wave of frequency (wq -|- traversing the medium. 



/35a;=/3\/l 



(112) 



The power spectrum of radiation is obtained from a{6ui) 
by taking its squared modulus. 

The mode sum of Eq. (Ill) may be evaluated numer- 
ically. However, because of its resonant character which 
implies significant contribution from only those modes 
that have propagation constant close to the plane-wave 
propagation constant Ps^, accurate analytical approxi- 
mations are also possible. First of all, for large spheres, 
^ >> 1, since all of the mode propagation constants 
7„ are largely real with relatively small imaginary parts, 
there is little radiation at frequency detunings that lie in 
the range (0, r?.o|/ip/(3fi-)) for which Psu) is purely imagi- 
nary. Due to a strong resonant interaction of light with 
matter, radiation simply cannot propagate far in this fre- 
quency range. There is thus a gap in the spectrum in this 
range, a feature of resonant radiation that was already 
noted in the one-dimensional problem [1]. 

Outside this gap, the propagation constant jSg^ is ei- 
ther large or small compared to the radius parameter 
and the following excellent analytical approximations can 
be verified for the propagation constants 7„ of the signif- 
icantly contributing modes, those for which Psw- 



=lnR 

(n-i)7r 

.riTT 
n-K — I — 



for 5lj <0 



for Sco > no|/ip/(3/i). 



(113) 



When these expressions are substituted in the resonant 
denominators of the terms in the mode sum (111) and in 
the large- argument approximation for ji(7„i?), namely 
— cos(7„i?)/(7„i?), and 7„i? is replaced by 13sm in terms 
that vary slowly from one mode to the next, the following 
analytical forms result for the spectral amplitude: 
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X < 



(-1)" 



y (-1)" 



for Suj <0 

for 5w > no|/Ip/(3fi), 



(114) 



where / = ro/i? denotes the fraction of the spherical ra- 
dius initially excited and 75;^ an effective propagation 
constant that takes slightly different forms below and 

above the gap, 



7. =^ 



for 5uj < 



■i/P) for 5u; > no 1/217 (Sfi.). 



(115) 



An overall factor that differs from 1 by terms of order 
1/(3 has been dropped from Eq. (114) which is bound 
to be accurate for large values of /3 only. Both sums in 
Eq. (114) may be evaluated in closed form. The second 
sum is well known [18] . The evaluation of the first sum by 
means of contour integration is presented in Appendix B. 
The following closed-form expressions are then obtained 
for the spectral amplitude below and above the frequency 
gap: 



a{Scj) 



h.2 2 

R 



X < 



1/3 

1 



3 Suj/Ssu 



1 



cos JSuj 

1 1 



for 5uj<0 



75a 



for 6u! > 



3h 



(116) 



The squared modulus S{5oj) = |a(5ti;)p represents the 
power spectral density of radiation. 

We plot in Figs. 4 the power spectrum S{5uj) for a 
variety of values of /3 and the fractional excited radius 
/ = ro/R- For small spherical samples with linear di- 
mensions comparable to the wavelength of light, as for 



/3 = 1, coherent radiation initially proceeds by the fast 
decaying superradiant mode we discussed in Sec. II, but 
the relatively slowly decaying modes that we discussed in 
Sec. V.B.I continue to radiate for long periods of time. 
Because of the simple proportionality of the decay rate 
to the width of the spectrum radiated by a mode, the su- 
perradiant mode furnishes a broad spectral background 
on which are coherently superposed sharper line spec- 
tra corresponding to the slowly decaying modes. Each 
spectral peak corresponding to a mode is centered at a 
frequency detuning equal to the imaginary part of the 
decay constant of that mode. We display the power spec- 
trum of emitted radiation for four different values of the 
fractional excited core radius, /, namely 0.1, 0.5, 0.75, 
and 1, in Figs. 4(a)-(d). Since the more localized ini- 
tial excitations are made up of a larger number of weakly 
decaying modes, all superposed coherently with the fun- 
damental supperradiant mode, the overall power spec- 
trum consists of a broad peak and a fine structure of 
ever narrower peaks that accumulate below zero detun- 
ing, as seen most dramatically in Fig. 4(d). An initial 
excitation of the full sphere corresponds, by contrast, to 
only a small admixture of the weakly decaying modes, 
that are visible in Fig. 4(a) as small peaks on the broad 
background provided by the superradiant mode. When 
considered in sequence, these four subplots also illustrate 
how even for small values of /3 a frequency gap develops 
over the interval [0, no/x^/(3fi)] as the initial excitation 
of the sphere is more and more tightly confined close to 
its center. The broad background contribution of the su- 
perradiant mode gets progressively smaller, yielding little 
power at any positive frequency detuning including this 
gap. 



For larger spheres too, a similar fine structure is ob- 
tained in the spectrum, but the emergence of a well de- 
fined frequency gap we discussed earlier is unmistakable 
with j3 increasing in value from 10 to 100, as we see in 
Figs. 5 and 6. The more localized the initial excitation 
the more prominent the peaks to the left of the gap. This 
is a result of the fact that the peaks with negative detun- 
ing correspond to modes that have propagation constants 
that are large compared to ko and which are therefore 



needed to make up an excitation that is localized on the 
sub- wavelength scale corresponding to /3/ << 1. The an- 
alytical approximation (116), shown by a dashed curve 
on each plot, is already accurate for /3 = 10, and is nearly 
indistinguishable from the numerically exact results for 
/3 = 100. For /3 = 10 the most superradiant mode has a 
decay rate that, according to expression (78), is a factor 
ln/3 = 2.3 times smaller than nojj? (3 / {2>h). This accounts 
for the considerably narrower background provided by 
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this mode to the emission spectrum in Fig. 5(a) than 
that present in Fig. 4(a) for (3 = 1, even though that 
mode dominates the other modes when the spliere is uni- 
formly excited. For much larger values of /3, as in Figs. 6, 
the superradiant modes, which are of order (ln/3)/7r in 
number, each correspond to a spatially nonuniform ex- 



citation given by Eq. (53) with ^ p. A uniformly 
excited core, regardless of its fractional radius, /, is thus 
comprised of large numbers of superradiant and weakly 
decaying modes, which always yield a rich spectrum of 
narrow peaks that accumulate on either side of the fre- 
quency gap. 



B. Time Dependence of the Radiated Power 

The power dW/dt radiated by the resonantly excited 
sphere is given by integrating the normal component of 
the Poynting vector over the surface of the sphere. For 
/3 >> 1, this integral may be reduced to a rather simple 
final form, as we show in Appendix C, 



dW 
~dt' 



kl{l + P') 



\d{R,t)\' 



(117) 



where d(R, t) is the amplitude function that determines 
the electric displacement field at the spherical surface, 
r — R, via relation (104). 



1. Radiation from a Large Sphere, j3 » \, with a Small 
Excited Gore, / << 1 

The early stages of radiation are dominated by the fast- 
decaying, superradiant modes, while the excitation resid- 
ing in the more weakly radiating modes is slow to escape 
the medium. The latter are accompanied by an oscilla- 
tory exchange of energy between the atoms and the ra- 
diation field, as we shall see presently. Similar attributes 
of cooperative emission from extended media have also 
been seen for a symmetric Dicke one-photon excitation 
[14, 16]. 

If we use in expression (106) the approximate eigen- 
value formulas (77) and (78), valid for the superradiant 
modes of a large sphere, and then evaluate the sum (104), 
we arc assured of good accuracy for the times when those 
modes are the ones actively radiating. Since the super- 
radiant modes have propagation constants that do not 
differ much from the free-space one, fco, we first approxi- 
mate (Ajq — 7^) by 2fco(fco — 7„), 7^ by k^, and then substi- 
tute formulas (77) and (78) in the expression (106). We 
assume here a small initially excited core, ro << R, so we 
may replace ji (7n?'o) by ji {Pf), while the large- argument 
approximation for ji (7n-R) coupled with expression (77) 
reduces it to the form 



MlnR) 



COS JnR 

InR 



(-1) 



n+1 ' 



t/4 



VP 



(118) 



where we used /3 >> 1 to ignore a term of relative or- 
der 1/p and replaced 7„i? in the denominator simply by 
p. With these approximations and upon extending the 
sum (104) to also include all negative integral values of 



the index n, which adds spurious terms to the sum that 
we later subtract out approximately, we may express the 
field amplitude d{R, t) at the surface of the sphere as 



d{R, t) 



E 



fp'hiPm 

(^_l^n^iat/[nn — P — (i/2) ^^P] 



[nn- P - {i/2) In P] 



(119) 



This sum may be evaluated exactly by the method of 
contour integrals in the complex plane, as shown in Ap- 
pendix D, and the following asymptotic expansion in 
powers of 1 //3 is obtained: 

d{R,t) =2e'-/^^fp'MPf)Po^ 

^E( — ) M'Wi'2n + l)at), (120) 



n=0 



where Jo denotes Bessel function of the first kind of order 

0. For large values of /3, the first term alone suffices to 
furnish an accurate result for d{R, t) and thus for the 
radiated power (111), 



^ « -^c.ofcor^[ii(/3/)]^|Po|Vo^(2v^). 



(121) 



The oscillatory time dependence (121) of the radiated 
power represents an exchange of energy between the field 
and the polarization of the radiating medium. Such os- 
cillatory energy exchange is characteristic of any radia- 
tion problem in which many modes of comparable de- 
cay rates that are detuned by different amounts from the 
resonance frequency participate coherently. A squared 
Bessel function time dependence related to Eq. (121) 
was first derived by Burnham and Chiao [19] in the 
one-dimensional context of radiation from a semi-infinite 
medium that is coherently generated in the wake of a 
sweeping (5— function excitation pulse. 

An expression for the radiated power that is somewhat 
more accurate than Eq. (121) is given by subtracting out 
the spurious terms, those with n running from to —00, 
that we included in the sum (119) in order to derive 
the Bessel- function result (120). Because of their non- 
resonant character, the sum of these spurious terms can 
be evaluated approximately, as we show in Appendix E, 
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and subtracted from Eq. (119). This procedure yields the 
following two-term result for the radiated power: 



dW _ IGtt 



at 



Jo(2\/^)sin(/3 + at//3)e 



(122) 



In Figs. 7-9, wc plot the power radiated by a spherical 
medium for the same values of the radius parameter /3 for 
which the radiated spectrum was considered in Figs. 4-6. 



For a uniformly excited sphere (/ = 1) with (3 = 1, the 
fundamental superradiant mode is nearly the only one 
excited, which implies a purely exponential decay of the 
radiated power, as seen in Fig. 7(a). When the medium 
is initially excited from the center out to only a tenth of 
the radius, as in Fig. 7(b), a significant number of weakly 
decaying modes of comparable magnitude are present as 
well. The initial precipitous drop in power results from 
the decay of the superradiant mode, but the subsequent 
emission has an oscillatory time dependence due to radi- 
ation emitted coherently by the weakly decaying modes. 
The oscillations are slow with a long quasi-period, due 
to rather small differences in the frequency detunings of 
these modes, as we noted earlier. 



As we increase the value of (3, we find inevitably the 
emergence of an oscillatory time dependence for power 
that for early times is described approximately by the 
Bcssel hmction result (122). This is particularly accu- 
rate for the larger value, 100, of /3 when the fractional 
excited core radius, /, is small, as seen in Figs. 9(b). 
The approximate result (122) ceases, however, to agree 
with the numerically exact behavior at long times due 
to the fact that the former wrongly presupposes that the 
superradiant modes continue to dominate the radiation 
for all times. In reality, the weakly decaying modes con- 



tinue to radiate for long times, well after the superradiant 
modes have radiated away nearly all of their excitation. 
The differences between the frequencies of neighboring 
modes, that determine the detailed character of the tem- 
poral oscillations, are dissimilar for the two varieties of 
modes, accounting for the failure of the approximate re- 
sult (122) to describe the long-time behavior of the radi- 
ated power. It is worth noting that times greatly exceed- 
ing those plotted would be necessary to see the stretched 
exponential decay of polarization energy predicted in the 
one-dimensional problem of localized excitation [1] . 



VIII. CONCLUDING REMARKS 

Coherent radiation from a sphere of polarizable atoms 
with a single resonant excitation energy level can exhibit 
a rich variety of spectral and temporal characteristics. A 
fully excited sphere with uniform polarization radiates 
much like a single atom, when its radius is small com- 
pared to the wavelength of emission. When coherently 
excited, a large sphere, by contrast, almost always radi- 
ates a superposition of strongly decaying, or superradi- 
ant, and weakly decaying, or subradiant, modes. These 
exponentially decaying modes may be classified as mag- 
netic and electric multipole modes of various angular mo- 
mentum orders. A particularly simple, yet sufficiently 
general, case of our problem is that of radiation from 
a uniformly excited concentric spherical core within an 
otherwise unexcited spherical medium. We have demon- 
strated that such an excitation radiates as a pure electric 
dipole, regardless of the size of the medium or its excited 
core. This specific case has been studied in detail in the 
present paper. 

The characteristic differences between the frequency 
detunings and decay rates of the modes give rise to ra- 
diation that has a wealth of sharp peaks and valleys in 
its spectrum. It also has a frequency range in which 



radiation cannot propagate far because of its strong res- 
onant interaction within the medium, and is thus un- 
able to escape it. We derived an analytical approxima- 
tion of this spectrum that involves simple trigonometric 
functions and is highly accurate for large values of /?. 
The time dependence is correspondingly quite involved, 
with a complicated oscillatory behavior that persists for 
long times. For times that are not too long, a simple 
analytical result based on the inclusion of superradiant 
modes alone furnishes a good approximation to the ex- 
act time dependence. Many of these characteristics of 
emission, previously noted in our treatment of the one- 
dimensional coherent radiation problem as well, are likely 
to survive the change of geometry, provided geometrical 
length scales remain large compared to the characteristic 
wavelengths of emission. A well localized excitation deep 
in the interior of an extended medium of arbitrary geom- 
etry will, for example, remain trapped for long periods 
of time, releasing energy only slowly, unless incoherent 
processes intervene. 
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Appendix A: Evaluation of the Inner Products in 
Eq. (99) 

When expression (102) for d{r,0) is substituted into 
the inner product formula (86), we obtain 

(d(r,0),ji(7„r)) = -4fc3r2i 

X {Rji{koro)h['\koR)[Rji{jnR)]' 

+ ll f r^h^i^ {kor> )ji {kor< )h Mdr] , 

(Al) 

where a prime superscript denotes a derivative with re- 
spect to the radial coordinate, here R. Since r> is the 
larger of r, tq and the smaller, to evaluate the integral 
in Eq. (Al), we write it as a sum of two integrals 

hiikoro) / r^ji{kQr)ji{jnr)dr 



Bessel differential equations, 

[rze{ar)]" + [a^ - £{£ + l)/r'^][rze{ar)] = 0, 



(A3) 



[rwiipr)]" + [p^ - £{£ + l)/r'^][rwi{pr)] = 0. (A4) 

By multiplying Eq. (A3) by rwe{pr) and Eq. (A4) by 

rz£{ar), then subtracting one resulting equation from the 
other, and finally integrating both sides over r followed 
by a simple rearrangement of terms, we obtain 



/ 



Zi{ar)'we{Pr)r^ dr = 



(a2 - /32) 

X J dr[rze{rwe)" — rwe{rzi)"]. 

(A5) 



Since the integrand of the right hand side of Eq. (A5) 
is the derivative of [rz^(ru>^)' — rw^(rz^)'], its integral is 
trivially evaluated, and Eq. (A5) reduces to the form 



+ ji{koro)[ r'^h^i\kor^ji{'ynr)dr. (A2) / r'^ ze{ar)we{l3r)r'^ dr = ^ \TZi{rwi)' -rwi{rzi)'\ 

Jrn J (a - P ) 



Indefinite integrals of the form J Zi{ar)wi{j3r)r'^dr, 
where z^ and are any two spherical Bessel functions 
of order I, can be computed analytically. To see how, 
we first recognize that Bessel functions obey appropriate 



(A6) 

The indefinite-integral formula (A6) may now be used 
to compute the two definite integrals in expression (A2), 
and thus the integral in Eq. (Al) for which the following 
expression results: 



{ji(fcoro) [Rh^^\koR){Rh{^nR))' - Rh{lnR){Rh^P{koR))' 



jiilnTo) roji{koro)iroh\^\koro)y - ro/ii^^(fcoro)(roji(fcoro))'j |. 



(A7) 



Use of a Wronskian identity turns the terms within the 

second pair of brackets in Eq. (A7) into i/fco, while the 
terms within the first pair of brackets can be combined in 
view of the eigenvalue relation (64). These simplifications 
reduce Eq. (A7) to the form 



r'^h^i\kor>)ji{kor<)ji{'ynr)dr = 



1 



X {ji{koro)Rh['\koR)[Rn{7nR)]'{l - kl/jl) 

{i/ko)Mlnro)}- (A8) 



+ 



When the integral (A8) is substituted into the right hand 
side of Eq. (Al), the inner product attains its final form 



(d(r,0),ji(7„r)) = -4ikyo\ —Po 



3 ""{kl-iD 



jiikoTo). 
(A9) 



The inner product {ji{'yn'r'),ji{'ynf)) involves, as defi- 
nition (86) indicates, the integral 



drjl{'^nr)dr, 



which can be evaluated exactly in closed form by means 
of the indefinite-integral identity [17] 



/ 



x'^jl{x)dx = Y[ii(a;) - jo(x)j|(x)]. 



This yields the following exact result for the inner prod- 
uct: 



_2 r 

(ii(7nr-),ii(7nr-)) = — \xji{x)[xji{x)]' 



y[jtW-Jo'(^)j2(^)]}, 

(AlO) 
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where x = "fnR- By substituting the expUcit trigono- 
metric forms for the spherical Bessel functions of orders 
0,1,2, 



jo{x) = , ji{x) 



sm X cos X 



j2ix) 



n 3 

sm a; ^ cos x, 



x-' X 



and performing simple algebraic manipulations, we may 
reduce Eq. (AlO) to the form 



(il(7n7-),jl(7n0) = 



-1 



a; — I 1 ^ I sm X cos x 

x-' 



cos X 



(All) 



For large \x\ ~ |7,i|-R, only the first term in the square 
brackets in Eq. (All) is important, 

{jiMJiM) « -R, bn\R » 1, (A12) 

which is the expression used in arriving at Eq. (106). 
I 



Appendix B: Evaluation of the First Sum in 
Eq. (114) 



Consider the contour integral 



I{a) 



dz 



Q z{z — a) cos^' 



(Bl) 



where the contour C may be taken to be a circle of ra- 
dius Ntt in the complex-z plane. Take A'^ to be a positive 
integer and a to be a finite complex number. In the limit 
N — >■ oo, the integral (Bl) must vanish, since the inte- 
grand goes to zero faster than 1 /N'^ while dz grows only 
linearly with A^. But, by the residue theorem, the inte- 
gral on the RHS is simply 27ri times the sum of residues 
of the integrand at all of its poles in the finite complex 
plane. In the limit A — > oo, the integrand has only sim- 
ple poles at 0, a, and (n— l/2)7r, n = 0, ±1, ±2, . . ., where 
the residues are easily evaluated, and the following sum 
formula results: 



= 27r J — 



a a cos a 



E 



1 



(n - l/2)7r[(n - l/2)7r - a] sm{n - l/2)7r j " 
I 



(B2) 



Note that sin(n — l/2)7r = (—1)" ^. The infinite sum by (1 — n') in the part of the sum that is over n = — oo to 
may be re-expressed as a one-sided sum by relabeling n and then dropping the prime from n', as shown below: 



E 



(-ir 



(-ir 



{n-l/2)[{n-l/2)Tr- 



E (n_i/2)[(n-l/2)7r. 



E 



(-1) 



l-n' 



^^(l/2-n')[(l/2-n')7r- 



(-ir 



6l(--V2) 



2a y: 



(n-l/2)7r-a {n-l/2)7T + a 
(-1)" 



(n-l/2)[(n- 1/2)2^2 _ ^2]- 
I 



When this result is substituted into Eq. (B2), we obtain 
a closed-form expression for the sum needed in Eq. (114), 
namely 



Appendix C: Cycle- Averaged Radiated Power 

When averaged over the fundamental oscillation pe- 
riod, 27r/aJoj the Poynting vector takes the form 



(-1)" 



S{r,t) = 2c Re E{r,t) x B*{f,t) 



(CI) 



1 - 



„_i ~ l/2)'''[(^ ~" l/2)^7r^ — a^] 2a^ \ cos a J integral of the normal component of S over the sur- 

(B3) face of the sphere, r = R, gives the total cycle-averaged 
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power, —dW/dt, radiated by the sphere at time t, 

-^=Rjd^nf-S\r=H. 



(C2) 



The electric field E is given by the expression (107)), 
while the magnetic field B, which obeys the Maxwell 
equation 



D = -^V X B, 

ko 



may be read oB from Eq. (103) for D, 

1 



B{r,t) 



2ik, 



■d{r,t)LYw. 



(C3) 



Using Eq. (CI) and a simple vector triple product rear- 
rangement, we may write r • 5 as 

f-S = 2cRe{fxE)-B*. (C4) 

In view of the form (107) for E, the vector identity 

d 

fx {V X A) = V(f • A) - A-r—A, 

and the operator identity r • Z = 0, it follows that 

^xE = -^^§^[re{r,t)]LY,o{ii). (C5) 

With the help of Eqs. (C5) and (C3) and noting that 
L* = —L, we may reduce Eq. (C4) to the form 



r-S = Re 



7C -» -» -* 

^d*{r,t)[re{r,t)]'LYw{n) ■ LY^m 



(C6) 

Integrating Eq. (C6) over all solid angles with the help 
of a normalization integral for spherical harmonics, 



j d^n LYio ■ LYio = -2, 



produces the following result for cycle-averaged power 
(C2) at the spherical surface, r = R: 



dW 
'~dt' 



-2i?Re 



tc 

—^^d*{R,t)[R e{R,t)]' 



(C7) 



If we now use expansions (104) and (108) and make use 
of the first equality in Eq. (109), we may write Eq. (C7) 
as the double mode sum 



(6-^**6-^'*}. (C8) 



Use of the eigenvalue equation (64) and the easily derived 
equality 



1, 



helps us express the power —dW/dt in the simple form 
-Itt = fc4(rT^lE^"Ji(7"^)exp(-A„i)| 

= '-^\d{r,t)\\ (C9) 
valid when ^ >> 1. 



Appendix D: Evaluation of the Sum in Eq. (119) 

Sums of this form are most simply evaluated by inte- 
grating the associated function 



^ ^iat/(-nz-u) 

sin TTZ {nz — u) ' 



(Dl) 



where u = /3+ {i/2) In /3, over a closed contour at infinity 
that encloses all of the simple poles at all integers, z = n, 
n = 0,±1,±2,..., and the essential singularity at 2; = 
u/tt of the function f{z). The integral over the contour 
vanishes, since the integrand f{z) decays to sufSciently 
rapidly as \z\ — > 00. From the residue theorem of analytic 
function theory, then, the sum of the residues at the poles 
and at the essential singularity must also vanish. Since 
the residue at the pole 2; = n is just 



(nTT — u) 

the following sum formula results: 



E 



^iat / {mr—u) 

{rnr — u) 



Residue 



^iat/ {z — u) 



sin 770 {z — u) 

(D2) 

The residue at the essential singularity may be com- 
puted by expanding the exponential in f{z) in a power 
series and using the standard formula for the residue at 

a pole of arbitrary order in each power-series term. This 
procedure leads to the following sum formula: 



E (-1)" 



(nTT — u) 



= -E 



n=0 



(n!)2 dit" 



smu 



(D3) 

The result (D3), although exact, is not particularly useful 
since it trades one sum for another. However, when ^ » 
1, it can be expanded in a series of terms which decrease 
in magnitude as increasing positive powers of 1//3. To do 
so, note that wc may write 



1 — 2iexp(m) 
sinu 1 — exp(2m) 



= —2i exp[i(2m -|- l)u]. 



m=0 
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and therefore 

d" 1 



du" sin u 



= -2i i"(2m + l)"exp[i(2m+l)M]. (D4) 



m=0 



Substitution of the result (D4) into Eq. (D3), followed by 
an interchange of the order of the m and n sums on the 
right-hand side (RHS), leads to the following asymptotic 
sum formula: 



- e-^/(— ") cxp[/(2m. + l)/3] 

^ ^ ' (nn-u) ^„ 

^ ' m=0 ' 

-(2m+ l)at]' 



n=0 



(D5) 



where the value, cxp(±iu) = cxp(±i/3)(/3)T^/^, was used 
to replace u in terms of ,5. Since the n-sum on the RHS 
of Eq. (D5) is simply ;i pcn\ cr-scries expansion of the 
Bessel function Jo(2-\/ (2m + l)at), Eq. (D5) reduces to 
a simpler form, 



(nTT — w) 



2^E 



m=0 



exp[i(2TO + l)/3] 

^m+l/2 



X Jo(2v'(2m+ l)ai). 



(D6) 



where the sum is now only over even non-positive inte- 
gers n: n — 2m, m = 0,-1,-2,... Since the sum on 
the right hand side of Eq. (El) may be regarded, again 
approximately, as an integral over m, we may write it as 

5^ (-l)"5(n) « f dmg'{2m) = \ f dng'{n). 

n=-oo 

(E2) 

Since the integrand is a total derivative with respect 
to the integration variable n, its integral is trivial, and 
Eq. (E2) reduces to the simple form 



(E3) 



since g'(— oo) = as a necessary consequence of the con- 
vergence of the original sum. 

Use of this general result immediately proves the fol- 
lowing sum formula: 



E (-1)" 



^iat/ {wK — u) 

{rnr — u) 



^ g-iat/u 



(E4) 



Noting that u = 13 + (i/2)ln/3 and l/u « 1//3 - 
(i/2)(ln/3)//3^, valid for large /3, the preceding result is 
approximately the same as 



Appendix E: Approximate Evaluation of a Certain 
Sum 

When a function g{n) changes slowly from one integer 

value of n to the next, the sum X]n=-oc(^l)".9('^) be 
evaluated by noting that the sum of each successive pair 
of terms since they have opposite signs is approximately 
the same as the first derivative g'{n). We may thus write 



Y (-i)"5(«) ^ E ^'M' 



(El) 





E 

I — — oo 



(-ir- 



ciat/(nn-u) ^ g-iat/0g-atln^/{20^) 

(riTT - u) 2 ^ ■ 



(E5) 

Subtracting this result from the first term of Eq. (D6), 
which for large /3 gives the sum over all integer values 
of p for the same summand, taking the squared modulus 
of the resulting difference, and then keeping only the two 
most significant powers of l/\/^ yield the result (121) for 
the corresponding sum over only positive integral values 
of n. 
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FIG. 4. The frequency spectrum of the power radiated by 
the medium, in arbitrary units, as a function of the frequency 
detuning, in units of no |/ip/3/(3/i), for /3 = 1, with (a) / = 1 
(uniform initial excitation); (b) / = 0.75; (c) / = 0.5; and (d) 
/ = 0.1 
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FIG. 5. Same as Fig. 4(a), except (a) ^ = 10, / = 1 and (b) 
13 = 10, f = 0.01. 
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FIG. 6. Same as Figs. 5(a), except (a) P = 100, / = 0.01 and (b) /3 = 100, / = 0.0001. 




FIG. 7. The power radiated by the medium, in arbitrary 
units, as a function of time, in units of [no|/2|^y0/(3/i)]~^, for 
13 = 1 with (a) / = 1.0 and (b) / = 0.1. 




FIG. 8. Same as Figs. 7(a) and (b), except (a) /3 = 10, 
/ = 1.0 and (b) /3 = 10, / = 0.01. 




FIG. 9. Same as Figs. 7(a) and (b), except (a) ^ = 100, / = 0.01 and (b) /3 = 100, / = 0.001. 



